A New Metaheuristic Bat-Inspired Algorithm 

Xin-She Yang 
Department of Engineering, University of Cambridge, 
Trumpington Street, Cambridge CB2 IPZ, UK 



Abstract 

Metaheuristic algorithms such as particle swarm optimization, firefly algorithm and 
harmony search are now becoming powerful methods for solving many tough optimiza- 
tion problems. In this paper, we propose a new metaheuristic method, the Bat Algo- 
rithm, based on the echolocation behaviour of bats. We also intend to combine the 
advantages of existing algorithms into the new bat algorithm. After a detailed formu- 
lation and explanation of its implementation, we will then compare the proposed algo- 
rithm with other existing algorithms, including genetic algorithms and particle swarm 
optimization. Simulations show that the proposed algorithm seems much superior to 
other algorithms, and further studies are also discussed. 
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1 Introduction 

Metaheuristic algorithms such as particle swarm optimization and simulated annealing 
are now becoming powerful methods for solving many tough optimization problems [3- 
7,11]. The vast majority of heuristic and metaheuristic algorithms have been derived 
from the behaviour of biological systems and/or physical systems in nature. For exam- 
ple, particle swarm optimization was developed based on the swarm behaviour of birds 
and fish [71 [S] , while simulated annealing was based on the annealing process of metals 

New algorithms are also emerging recently, including harmony search and the firefly 
algorithm. The former was inspired by the improvising process of composing a piece 
of music [4] , while the latter was formulated based on the flashing behaviour of fireflies 
[14]. Each of these algorithms has certain advantages and disadvantages. For example, 
simulating annealing can almost guarantee to find the optimal solution if the cooling 
process is slow enough and the simulation is running long enough; however, the fine 
adjustment in parameters does affect the convergence rate of the optimization process. 
A natural question is whether it is possible to combine major advantages of these 
algorithms and try to develop a potentially better algorithm? This paper is such an 
attempt to address this issue. 
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In this paper, we intend to propose a new metaheuristic method, namely, the Bat 
Algorithm (BA), based on the echolocation behaviour of bats. The capability of echolo- 
cation of microbats is fascinating as these bats can find their prey and discriminate 
different types of insects even in complete darkness. We will first formulate the bat 
algorithm by idealizing the echolocation behaviour of bats. We then describe how it 
works and make comparison with other existing algorithms. Finally, we will discuss 
some implications for further studies. 

2 Echolocation of bats 

2.1 Behaviour of microbats 

Bats are fascinating animals. They are the only mammals with wings and they also have 
advanced capability of echolocation. It is estimated that there are about 996 different 
species which account for up to 20% of all mammal species [1] [2] . Their size ranges from 
the tiny bumblebee bat (of about 1.5 to 2g) to the giant bats with wingspan of about 2 
m and weight up to about 1 kg. Microbats typically have forearm length of about 2.2 to 
11cm. Most bats uses echolocation to a certain degree; among all the species, microbats 
are a famous example as microbats use echolocation extensively while megabats do not 

[mils]. 

Most microbats are insectivores. Microbats use a type of sonar, called, echolocation, 
to detect prey, avoid obstacles, and locate their roosting crevices in the dark. These bats 
emit a very loud sound pulse and listen for the echo that bounces back from the sur- 
rounding objects. Their pulses vary in properties and can be correlated with their hunt- 
ing strategies, depending on the species. Most bats use short, frequency-modulated sig- 
nals to sweep through about an octave, while others more often use constant-frequency 
signals for echolocation. Their signal bandwidth varies depends on the species, and 
often increased by using more harmonics. 

2.2 Acoustics of Echolocation 

Though each pulse only lasts a few thousandths of a second (up to about 8 to 10 ms), 
however, it has a constant frequency which is usually in the region of 25kHz to 150 
kHz. The typical range of frequencies for most bat species are in the region between 
25kHz and lOOkHz, though some species can emit higher frequencies up to 150 kHz. 
Each ultrasonic burst may last typically 5 to 20 ms, and microbats emit about 10 to 20 
such sound bursts every second. When hunting for prey, the rate of pulse emission can 
be sped up to about 200 pulses per second when they fly near their prey. Such short 
sound bursts imply the fantastic ability of the signal processing power of bats. In fact, 
studies shows the integration time of the bat ear is typically about 300 to 400 ^s. 

As the speed of sound in air is typically v — 340 m/s, the wavelength A of the 
ultrasonic sound bursts with a constant frequency / is given by 



which is in the range of 2mm to 14mni for the typical frequency range from 25kHz to 
150 kHz. Such wavelengths are in the same order of their prey sizes. 
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Amazingly, the emitted pulse could be as loud as 110 dB, and, fortunately, they 
are in the ultrasonic region. The loudness also varies from the loudest when searching 
for prey and to a quieter base when homing towards the prey. The travelling range of 
such short pulses are typically a few metres, depending on the actual frequencies P^ . 
Microbats can manage to avoid obstacles as small as thin human hairs. 

Studies show that microbats use the time delay from the emission and detection of 
the echo, the time difference between their two ears, and the loudness variations of the 
echoes to build up three dimensional scenario of the surrounding. They can detect the 
distance and orientation of the target, the type of prey, and even the moving speed of 
the prey such as small insects. Indeed, studies suggested that bats seem to be able to 
discriminate targets by the variations of the Doppler effect induced by the wing-flutter 
rates of the target insects [1] . 

Obviously, some bats have good eyesight, and most bats also have very sensitive 
smell sense. In reality, they will use all the senses as a combination to maximize the 
efficient detection of prey and smooth navigation. However, here we are only interested 
in the echolocation and the associated behaviour. 

Such echolocation behaviour of microbats can be formulated in such a way that it 
can be associated with the objective function to be optimized, and this make it possible 
to formulate new optimization algorithms. In the rest of this paper, we will first outline 
the basic formulation of the Bat Algorithm (BA) and then discuss the implementation 
and comparison in detail. 

3 Bat Algorithm 

If we idealize some of the echolocation characteristics of microbats, we can develop var- 
ious bat-inspired algorithms or bat algorithms. For simplicity, we now use the following 
approximate or idealized rules: 

1. All bats use echolocation to sense distance, and they also 'know' the difference 
between food/prey and background barriers in some magical way; 

2. Bats fly randomly with velocity at position with a fixed frequency /min, 
varying wavelength A and loudness Aq to search for prey. They can automatically 
adjust the wavelength (or frequency) of their emitted pulses and adjust the rate 
of pulse emission r e [0, 1], depending on the proximity of their target; 

3. Although the loudness can vary in many ways, we assume that the loudness varies 
from a large (positive) Aq to a minimum constant value Amin- 

Another obvious simplification is that no ray tracing is used in estimating the time 
delay and three dimensional topography. Though this might be a good feature for 
the application in computational geometry, however, we will not use this as it is more 
computationally extensive in multidimensional cases. 

In addition to these simplified assumptions, we also use the following approxima- 
tions, for simplicity. In general the frequency / in a range [/min , /max] corresponds to 
a range of wavelengths [Amin, Amax]- For example a frequency range of [20kHz, 500kHz] 
corresponds to a range of wavelengths from 0.7mm to 17mm. 

For a given problem, we can also use any wavelength for the ease of implementation. 
In the actual implementation, we can adjust the range by adjusting the wavelengths 
(or frequencies), and the detectable range (or the largest wavelength) should be chosen 
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Bat Algorithm 

Objective function /(x), x = (.ti, XdY 
Initialize the bat population Xj (z = 1, 2, n) and 
Define pulse frequency fi at Xi 
Initialize pulse rates ri and the loudness Ai 
while (t <Max number of iterations) 
Generate new solutions by adjusting frequency, 

and updating velocities and locations /solutions [equations (0j to ^] 
if (rand > ri) 

Select a solution among the best solutions 

Generate a local solution around the selected best solution 

end if 

Generate a new solution by flying randomly 
if (rand < Ai k. /(x^) < /(x*)) 
Accept the new solutions 
Increase ri and reduce Ai 
end if 

Rank the bats and find the current best x* 
end while 

Postprocess results and visualization 



Figure 1: Pseudo code of the bat algorithm (BA). 



such that it is comparable to the size of the domain of interest, and then toning down 
to smaller ranges. Furthermore, we do not necessarily have to use the wavelengths 
themselves, instead, we can also vary the frequency while fixing the wavelength A. This 
is because A and / are related due to the fact A/ is constant. We will use this later 
approach in our implementation. 

For simplicity, we can assume / e [0,/max]- We know that higher frequencies have 
short wavelengths and travel a shorter distance. For bats, the typical ranges are a few 
metres. The rate of pulse can simply be in the range of [0, 1] where means no pulses 
at all, and 1 means the maximum rate of pulse emission. 

Based on these approximations and idealization, the basic steps of the Bat Algorithm 
(BA) can be summarized as the pseudo code shown in Fig. [TJ 



3.1 Movement of Virtual Bats 

In simulations, we use virtual bats naturally. We have to define the rules how their 
positions Xj and velocities V; in a d-dimensional search space are updated. The new 
solutions X* and velocities v* at time step t are given by 

fi — /min (/max /min)/^i (2) 



+ (x*-x,)/„ (3) 



(4) 

where (3 G [0, 1] is a random vector drawn from a uniform distribution. Here x* is the 
current global best location (solution) which is located after comparing all the solutions 
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among all the n bats. As the product Xifi is the velocity increment, we can use either fi 
(or Xi ) to adjust the velocity change while fixing the other factor Xt (or fi), depending 
on the type of the problem of interest. In our implementation, we will use /min = and 
/max = 100, depending the domain size of the problem of interest. Initially, each bat is 
randomly assigned a frequency which is drawn uniformly from [/min, /max]- 

For the local search part, once a solution is selected among the current best solutions, 
a new solution for each bat is generated locally using random walk 

Xncw = Xold + e^*, (5) 

where e € [—1, 1] is a random number, while A* —<Al > is the average loudness of all 
the bats at this time step. 

The update of the velocities and positions of bats have some similarity to the pro- 
cedure in the standard particle swarm optimization as fi essentially controls the 
pace and range of the movement of the swarming particles. To a degree, BA can be 
considered as a balanced combination of the standard particle swarm optimization and 
the intensive local search controlled by the loudness and pulse rate. 

3.2 Loudness and Pulse Emission 

Furthermore, the loudness Ai and the rate of pulse emission have to be updated 
accordingly as the iterations proceed. As the loudness usually decreases once a bat has 
found its prey, while the rate of pulse emission increases, the loudness can be chosen 
as any value of convenience. For example, we can use Aq — 100 and ^min — 1. For 
simplicity, we can also use Aq = 1 and Amin = 0, assuming ^min — means that a bat 
has just found the prey and temporarily stop emitting any sound. Now we have 

Al+'=aAl r*+i=r°[l-exp(-7t)], (6) 

where a and 7 are constants. In fact, a is similar to the cooling factor of a cooling 
schedule in the simulated annealing For any < a < 1 and 7 > 0, we have 

Al 0, r* -> r°, as t 00. (7) 

In the simplicity case, we can use a — "f, and we have used a = 7 = 0.9 in our 
simulations. The choice of parameters requires some experimenting. Initially, each 
bat should have different values of loudness and pulse emission rate, and this can be 
achieved by randomization. For example, the initial loudness Ai can typically be [1,2], 
while the initial emission rate can be around zero, or any value G [0, 1] if using 
([6]). Their loudness and emission rates will be updated only if the new solutions are 
improved, which means that these bats are moving towards the optimal solution. 

4 Validation and Comparison 

From the pseudo code, it is relatively straightforward to implement the Bat Algorithm 
in any programming language. For the ease of visualization, we have implemented it 
using Matlab for various test functions. 
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Figure 2: The paths of 25 virtual bats during 20 consecutive iterations. They converge into 
(1,1)- 



4.1 Benchmark Functions 

There are many standard test functions for validating new algorithms. In the current 
benchmark validation, we have chosen the well-known Rosenbrock's function 



/(x) = ^(1 - xff + 100(a;i+i - x^)^, -2.048 < a;, < 2.048, 



(8) 



and the eggcrate function 

5(2;, y) = + + 25(sin^a; + sin^ y), (a;, y) e [-27r, 27r] x [-27r, 27r]. (9) 

We know that /(x) has a global minimum fmm = at (1, 1) in 2D, while g{x,y) has a 
global minimum gmin = at (0, 0). De Jong's standard sphere function 



/i(x)=^x2, -10<a;, <10, 



(10) 



has also been used. Its minimum is ft-min = at (0, 0, 0) for any d > 3. 

In addition, we have also used other standard test functions for numerical global 
optimization fTU| such as Ackley's function 



s(x) = 20 + e-20exp 



0.2, 



1 1 

\ i=l i=l 
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Figure 3: The eggcrate function (left) and the locations of 40 bats in the last ten iterations 
(right). 



where —30 < xi < 30. It has the global minimum s„ 
Michalewicz's test function 



at (0,0,..., 0). 



/(x) = - XI ^"^(^ 



sin(^) , (m=10), 

TT 



(12) 



has d\ local optima in the the domain < a;i < tt where i = l,2,...,(i. The global 
minimum is /* « —1.801 for d — 2, while /, « —4.6877 for d ~ 5. 

In our implementation, we use n = 25 to 50 virtual bats, and a = 0.9. For Rosen- 
brock's 2-D banana function, the paths of 25 virtual bats during the consecutive 20 time 
steps are shown in Fig. [5] where we can see that the bats converge at the global opti- 
mum (1, 1). For the multimodal eggcrate function, a snapshot of the last 10 iterations 
is shown in Fig. |31 Again, all bats move towards the global best (0,0). 



4.2 Comparison with Other Algorithms 

In order to compare the performance of the new algorithm, we have tested it against 
other heuristic algorithms, including genetic algorithms (GA) [5l[6j[Tl], and particle 
swarm optimization (PSO) [7l[8]. There are many variants of PSO, and some variants 
such as the mean PSO could perform better than the standard PSO f^; however, the 
standard PSO is by far the most popularly used. Therefore, we will also use the standard 
PSO in our comparison. 

There are many ways to carry out the comparison of algorithm performance, and 
two obvious approaches are: to compare the numbers of function evaluations for a 
given tolerance or accuracy, or to compare their accuracies for a fixed number of func- 
tion evaluations. Here we will use the first approach. In our simulations, we use a 
fixed tolerance e < 10~^, and we run each algorithm for 100 times so that we can do 
meaningful statistical analysis. 
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For genetic algorithms, we have used the standard version with no ehtism with 
the mutation probabihty of Pm — 0.05 and crossover probabihty of 0.95. For particle 
swarm optimization, we have also used the standard version with learning parameters 
a — =2 and the inertia function / = 1 [71 [S]. The simulations have been carried out 
using Matlab on a standard 3GHz desktop computer. Each run with about 10,000 
function evaluations typically takes less than 5 seconds. Furthermore, we have tried to 
use different population sizes from n = 10 to 250, and we found that for most problems, 
n = 15 to 50 is sufficient. Therefore, we use a fixed population n = 40 for all simulations. 
Table 1 shows the number of function evaluations in the form of mean ± the standard 
deviation (success rate of the algorithm in finding the global optima). 

Table 1: Comparison of BA with GA, and PSO. 
Functions /Algorithms GA PSO BA 



Multiple peaks 
Michalewicz's (d=16) 
Rosenbrock's ((i=16) 
De Jong's (d=256) 
Schwefel's ((i=128) 
Ackley's (d=128) 
Rastrigin's 

Easom's 
Griewangk's 
Shubert's (18 minima) 



52124 ± 3277(98%) 
89325 ± 7914(95%) 
55723 ±8901(90%) 
25412 ± 1237(100%) 
227329 ± 7572(95%) 
32720 ± 3327(90%) 
110523 ± 5199(77%) 
19239 ± 3307(92%) 
70925 ± 7652(90%) 
54077 ± 4997(89%) 



3719 ± 205(97%) 
6922 ± 537(98%) 
32756 ± 5325(98%) 
17040 ± 1123(100%) 
14522 ± 1275(97%) 
23407 ± 4325(92%) 
79491 ± 3715(90%) 
17273 ± 2929(90%) 
55970 ± 4223(92%) 
23992 ± 3755(92%) 



1152 ± 245(100%) 
4752 ± 753(100%) 
7923 ± 3293(100%) 
5273 ± 490(100%) 
8929 ± 729(99%) 
6933 ± 2317(100%) 
12573 ± 3372(100%) 
7532 ± 1702(99%) 
9792 ± 4732(100%) 
11925 ±4049(100%) 



From Table 1, we can see that PSO performs much better than genetic algorithms, 
while the Bat Algorithm is much superior to other algorithms in terms of accuracy and 
efficiency. This is no surprising as the aim of developing the new algorithm was to try 
to use the advantages of existing algorithms and other interesting feature inspired by 
the fantastic behaviour of echolocation of microbats. 

If we replace the variations of the frequency fi by a random parameter and setting 
Ai = Q and = 1, the bat algorithm essentially becomes the standard Particle Swarm 
Optimization (PSO). Similarly, if we do not use the velocities, but we use fixed loudness 
and rate: Ai and r^. For example, Ai ~ Vi — 0.7, this algorithm is virtually reduced to 
a simple Harmony Search (HS) as the frequency /wavelength change is essentially the 
pitch adjustment, while the rate of pulse emission is similar to the harmonic acceptance 
rate (here with a twist) in the harmony search algorithm [4l I15| . The current studies 
implies that the proposed new algorithm is potentially more powerful and thus should 
be investigated further in many applications of engineering and industrial optimization 
problems. 

5 Discussions 

In this paper, we have successfully formulated a new Bat Algorithm for continuous 
constrained optimization problems. From the formulation of the Bat Algorithm and its 
implementation and comparison, we can see that it is a very promising algorithm. It is 
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potentially more powerful than particle swarm optimization and genetic algorithms as 

well as Harmony Search. The primary reason is that BA uses a good combination of 
major advantages of these algorithms in some way. Moreover, PSO and harmony search 
are the special cases of the Bat Algorithm under appropriate simplifications. 

In addition, the fine adjustment of the parameters a and 7 can affect the convergence 
rate of the bat algorithm. In fact, parameter a acts in a similar role as the cooling 
schedule in the simulated annealing. Though the implementation is more complicated 
than many other metahcuristic algorithms; however, it does show that it utilizes a 
balanced combination of the advantages of existing successful algorithms with innovative 
feature based on the echolocation behaviour of bats. New solutions are generated by 
adjusting frequencies, loudness and pulse emission rates, while the proposed solution is 
accepted or not depends on the quality of the solutions controlled or characterized by 
loudness and pulse rate which are in turn related to the closeness or the fitness of the 
locations/solution to the global optimal solution. 

The exciting results suggest that more studies will highly be needed to carry out 
the sensitivity analysis, to analyze the rate of algorithm convergence, and to improve 
the convergence rate even further. More extensive comparison studies with a more 
wide range of existing algorithms using much tough test functions in higher dimensions 
will pose more challenges to the algorithms, and thus such comparisons will potentially 
reveal the virtues and weakness of all the algorithms of interest. 

An interesting extension will be to use different schemes of wavelength or frequency 
variations instead of the current linear implementation. In addition, the rates of pulse 
emission and loudness can also be varied in a more sophisticated manner. Another 
extension for discrete problems is to use the time delay between pulse emission and 
the echo bounced back. For example, in the travelling salesman problem, the distance 
between two adjacent nodes/cities can easily be coded as time delay. As microbats use 
time difference between their two ears to obtain three-dimensional information, they can 
identify the type of prey and the velocity of a flying insect. Therefore, a further natural 
extension to the current bat algorithm would be to use the directional echolocation and 
Doppler effect, which may lead to even more interesting variants and new algorithms. 
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